o 
o 

j3 ; 

^ . Random Cluster Models on the Triangular Lattice 

O . 

L. Chayes and H. K. Lei 
February 2, 2008 

O 

S 

> 
> 



^Department of Mathematics, UCLA, Los Angeles, California, USA 



Abstract: We study percolation and the random cluster model on the triangular lattice with 
3-body interactions. Starting with percolation, we generalize the star-triangle transformation: 
C/3 , We introduce a new parameter (the 3-body term) and identify configurations on the triangles 

solely by their connectivity. In this new setup, necessary and sufficient conditions are found for 
I positive correlations and this is used to establish regions of percolation and non-percolation. 

, Next we apply this set of ideas to the q > 1 random cluster model: We derive duality rela- 

tions for the suitable random cluster measures, prove necessary and sufficient conditions for 
them to have positive correlations, and finally prove some rigorous theorems concerning phase 
transitions. 
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Introduction 



> 

CN . The study of duality relations for 2Z)-Potts systems is not a new topic. Indeed, it is older 

^ I than the model itself in the sense that [21] and [2] provided special cases long before the 

l/^ • general Potts spin-systems were formulated. While we will not dwell on the historical 

■ aspects of this subject, it is worth remarking that this line of study has had immeasurable 

impact on the entirety of two-dimensional statistical mechanics. Notwithstanding, the 
usual derivations of duality for Potts models (see [36] and references therein) suffer in three 
respects which we will describe in increasing order of importance: 
i • (I) There are informal aspects to many of the derivations and thus some effort - pre- 

Q I sumably small - would be needed to elevate these derivations to the status of mathematical 

O ' theorems. 

• (II) The various standard techniques, which include mapping to vertex models or the 
^ \ introduction of dual-spin variables in the form of constraints, do not include all relevant 

• values of parameters. In particular, the dual-constraints approach only makes sense for 



integer q > 2. It is only as an afterthought that duality relations for continuous g's are 
inferred from the analytic structure of the formulas produced for the integer q's. 

While we do not necessarily regard these two issues as being of great urgency, the third 
issue is considered to be pertinent both by mathematicians and physicists. 
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• (III) The result of a typical duality relation is the identification of the free energies at 
dual parameter values. Hence, as concerns the subject of phase transitions one is always 
left with an unsatisfactory provisional statement: //there is but a single non-analyticity, 
then this must occur at the self-dual point. 

It should be remarked that this third issue is certainly not "academic". In particular, 
in the so called rs-models [12], which are generalizations of the Ashkin-Teller and/or the 
g-state Potts models (with q = rx s) there is a self-dual line through an intermediate phase 
where, apparently, nothing of interest transpires; c.f. [32], [26] and [9]. 

As an alternative to the "usual methods" , it is possible to establish duality via graphical 
representations, in particular the FK-representation [15], whereby the duality shows up on 
the level of the representation itself. Duality in this context is akin to (and a generalization 
of) the elementary sort of duality found in Bernoulli percolation. Hence, using percolation 
based techniques, genuine irrefutable statements can be made concerning the presence of 
phase transitions at points of self-duality. For example, on the square lattice, duality of the 
random cluster models has been used to establish rather sharp theorems concerning their 
phase structures [3], [11]. 

In this work we will study the q'-state Potts models - and their associated random cluster 
representations - on the triangular lattice. For these problems, the derivation is considerably 
more intricate than the square lattice; one must first go through the intermediate honeycomb 
lattice. The inevitable consequence of this contortion is the production of extra correlations 
in the dual model. In the language of spin-systems, these correlations translate into the 
phrase "three body interactions" but we iterate that the phenomenon is quite general and 
occurs even for percolation {q = 1). Well known exceptions to this rule are (i) The Ising 
spin-system at all couplings and (ii) A special point, called the star-triangle point, where 
by a miracle, the correlations in the dual model vanish. Since the star-triangle point is also 
a point of self-duality, it is readily identified as the transition point as was done in [20]. 
However, to the authors' knowledge, it is only for the case of percolation ([33], [34]) that a 
rigorous theorem along these lines has been established. 

The perspective of this work is that since we are generically stuck with the additional 
correlations after duality, then they should be in the model from the outset. We find that 
with the additional freedom of "three body interactions" , duality becomes a straightforward 
map in a two-dimensional space that has a self-dual curve of fixed points. One of the points 
on this curve - and of no particular significance - is the star-triangle point. This general 
picture has been known (and under appreciated) for quite some time: Duality relations on 
the level of free energies are derived in [36] using the methods of [5] - here for integer q > 2. 
Additional results along these lines are obtained in [35] , [37] and [4] via relations to vertex 
models. A cornerstone of the former work is a graphical expansion akin to what is developed 
here. However, in these works the representation was only employed as an auxiliary device. 
The full potential for relating percolation phenomena in the graphical representations to 
phase transitions (as defined by other means) and the use of the interplay between direct 
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and dual representations to elucidate this phenomena was not exploited. 

Prom the perspective of rigorous analysis, a significant problem emerges at the outset. 
In particular, the sorts of additional correlations introduced are not necessarily positive 
correlations. E.g. for the spin-systems, the extra interactions are, as often as not, antifer- 
romagnetic. While this may or may not alter the nature of the transition, it is an enormous 
technical obstacle since nearly all probabilistic arguments concerning systems of this sort 
are based on the positivity of correlations. To overcome these difficulties, we must intro- 
duce a reduced state space for the graphical models wherein positive correlations can be 
re-established. Notwithstanding, our techniques do not cover the entirety of the self-dual 
curve but this could in principle be accomplished by an extension of our scheme. Further, to 
avoid technical complications we deal exclusively with the isotropic case whenever possible: 
A priori, all three edges of the triangle have the same probability of being occupied. One 
might also, with some effort, extend various results proved here to the anisotropic cases. 

The remainder of this paper is organized as follows: In section two, we examine the 
case of percolation where the necessity of introducing local correlations is underscored. 
Here the star-triangle duality is generalized and relatively complete results for the phase 
diagram are derived. In section three, we study this problem for the q > 1 random cluster 
models. The duality of [35] and [4] are derived by graphical methods and wc characterize 
the conditions for positive correlations. Finally, in section four, wc show that in the region 
where correlations are positive, there is a phase transition (or at least critical behavior) at 
all points of self-duality. 



2 Generalized Star— Triangle Relations — Percolation 
2.1 The Classical Star— Triangle Situation 

In order to motivate our work, we first briefly describe the classical star-triangle relation. As 
mentioned above, we will treat the isotropic case, so let p be the probability that a bond is 
occupied. Now on any given triangle there are eight possible configurations; we denote their 
respective probabilities by e (empty), s (single), d (double) and a (full). Thus, for example, 
s = p{l —p)'^. Under the usual sort of planar duality, the triangular lattice problem becomes 
a problem on the honeycomb lattice where we could also associate a bond probability e.g. 

= 1—p. Considering only connectivity properties and integrating OT.it the central vertex 
returns us to a problem on the triangular lattice (but with the triangles inverted) . Using e* , 
s*, d* and a* to denote probabilities of the corresponding configurations, we easily arrive 
at 

e* =p^ + ^p'^(l-p)^ (2.1) 

s*=p{l-pf, (2.2) \^ 

M*+a* = (l-pf. (2.3) [22^ 
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Ostensibly, one would like to define a p* such that the right hand sides of (2.1), (2.2), 
and (2.3) are, respectively, (1 -p*)"', p*(l — p*Y ^i^d (jp*)^ + 3(p*)^(l —p*). However, for 
general p G (0, 1), this cannot be done - there are just too many equations. Explicitly, if 
we try to force this sort of duality, this in turn forces p to a particular value which, in fact, 
is the one for which p = p*. To see this, if we substitute (2.1) into (2.2) we get, in the 
variables R = p/{l - p) and R* = p*/{l - p*), the equation RR*{R* + 3) = 1. But the 
similar procedure on (2.3) and (2.2) gets us RR*{R + 3) = 1 thence any non-trivial solution 
requires R = R*. At p* = p, we see that p must satisfy: 

/ - 3p + 1 = 0, (2.4) [23^ 

which is of course the self-dual point of the classic star-triangle relation. 

2.2 Introduction of Correlations 

Overall, the above situation is clearly not suitable for the development of a general theory 
of duality. Clearly, if we wish to salvage this situation, the next step would be to put 
in some sort of correlations. A manageable way to implement correlations - which has 
its analogs in physical systems, c.f. subsection 3.1 - is to introduce correlations within 
triangles but to keep separate triangles independent. (Here, of course, we refer only to "up- 
pointing" triangles; configurations on the "down-pointing" triangles will be determined from 
the former.) 

A secondary consequence of the above duality experiment (on a single triangle) is the 
observation that, when the rinse cycle is finished, the dual model does not really distinguish 
between the double and full configurations. This is due to the fact that all we track are 
connectivities between sites and, in both situations, the triangle is fully connected. 

In this spirit, we might as well confine all of our attention to the three types of con- 
figurations listed in (2.1), (2.2), and (2.3); e.g., we define our model to have only five 
configurations on each triangle, namely empty, three singles and a full. So (in the fully 
isotropic case) we have five parameters: e, s and a with e + 3s + a = 1. We state without 
proof the following proposition concerning this model on the triangular lattice: 

[T] Proposition 2.1. Consider the model on the triangular lattice in which configurations on 
the up-pointing triangles are independent and confined to empty, singles and full with re- 
spective probabilities e, s and a. Then this model is dual to the model with parameters e* , 
s* and a* which are given by 

e* =a (2.5) | 24a 

s* =s (2.6) [24b 

a* =e (2.7) [24^ 

In particular, the condition for self-duality is just a = e. 
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We make a simple observation which will be useful in the next subsection: 



Corollary 2.2. For the parameters a, e as above and for r G [0, 1], the curve a + e = r is 
invariant under the *-map. 

In order to translate all of this into a statement about percolation properties of the 
model we will need to establish some FKG-type properties of the system. Since separate 
up-pointing triangles are independent this amounts to a problem on a single triangle. Here, 
unfortunately, we must prove the result for the anisotropic case as it will be needed later. 
First we need some basic definitions. 

Definition 2.3. Let O be a probability space with probability measure P. Let A C ^ he 
an event and let a; G fi. Then the indicator function 1a is defined by 

1 otherwise. 

If / is a function on J7, then E(/), the expectation (or mean value) of / is defined to be 

E(/) = / f{co)dP{u). 

Finally, we say the functions / and g have positive correlations if 

Eifg) > E(/)E(5). 

2t Theorem 2.4. Consider the above described 5-state system realized as bond configurations 
on a triangle: Let [S]i, [Sja and [SJa denote the events that the three various sides of the 
triangle are the sole bonds occupied with [A] and [E] denoting the full and empty configura- 
tions. Let u denote a measure on this system and let us denote the respective probabilities 
of the above-mentioned by si, S2, S3, a and e. It is assumed without loss of generality that 
si > S2 > S3. Then the necessary and sufficient condition for v to have positive correlation 
is 

ae > si{s2 + S3) 

Proof: To prove the necessity of the condition ae > si{s2 + S3), note that if /(si) = 0, 
/(S2) = 1 = /(S3) = 1, /(a) = 1 and /(e) = 0, and g{si) = 0, g{s2) = gis^) = 1, g{a) = 1 
and g{e) = 0, then 'E{fg) > E(/)E(5) gives exactly that ae > si(s2 + S3). For sufficiency, 
we aim to show that 

E(/ff) - E(/)E(<7) > (2.8) ^ 

To simplify matters we first note that (2.8) is not changed by adding constants to / and 
g. Thus we may assume that / and g are overall non-negative and (by subtracting f{E) 
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and g{E) respectively) vanish on the lowest configuration. Similarly, the truth or falsehood 
of (2.8) is unaffected by the scaling of / and g so we may as well assume that /([A]) = 
g([A]) = 1. 

Next let (T be a permutation on three letters such that /([S]cti) >/([S]o-2) ^/([SJas)- Then 
we are down to six parameters: for convenience let Xi, X2, denote /([Sjo-j), /([S]ct2)) 
/([Sjo-j), respectively. Similarly define yi, y2, and 2/3 for g. We assume that some of these 
parameters are non-trivial, for otherwise the theorem is already proved. 

Next we observe that any increasing function is automatically positively correlated with 
1[A], the indicator of the top configuration. Indeed (with all of our simplifications enforced), 
E(l[/\](7) = a, whereas E(1^)E(5() = aE{g), which is smaller. Thus, the quantity E{fg) — 
E(/)E(5() will decrease if we subtract from / the function XIa with A > 0. However, in 
order to keep / increasing, the most we can subtract is A = 1 — max{xi, X2, X3} = 1 — xi, 
by assumption. Thus, after this subtraction and more rescaling, we have that xi = 1. 

Similar considerations show that min{xi, X2, x^} = 0. To sec this one first observe that 
g is always positively correlated with the function 1 — Then subtracting from / the 

function X3(l — 1[e]) (where by assumption X3 = min{xi,X2, X3}) and rescaling again gives 
the desired conclusion. 

Given all these simplifications, we now have E(/(jr) = a+Stj-^yi + Sa,^X2y2 and 'E{f )'E{g) = 
(o + So-i +So-2 3^2) (a + 2/1 + S(T2?/2 + 0-32/3 )• Since the goal is to show that ^{fg) > 'E{f)'E{g), 
we may assume that yi = and ys = 1, since the coefficient of yi in E{f)'E{g) is smaller 
than in 'E{fg) and ys does not even occur in 'Ei{fg). 

Next one can check that / is positively correlated with l[s]^2 ■'"[Sl'^3 '^i^Y' 
this observe that E((l[s]^2 + 1[S]<t3 ) ~ (^2'So-2 + '^)(S(T2 + '^0-3 + a + .Sq-^ + e) whereas 

E(/)E(l[s]^2 + 1[S]^3 + 1[A]) = {xsa2 +a+Sai){sa2+Sa3 +a), SO the difference is ae-s„^ {sa^ + 
Scts), which is positive by hypothesis. It is also easy to check that g — 2/2(15^2 '^^'^z '^^^ 
is still increasing. Also, note that if y2 was equal to one before the subtraction, then after 
the subtraction g = and again the conclusion of the theorem holds trivially, so we may 
as well assume 1/2 7^ 1- As before. Subtracting and renormalizing, we acquire y2 = 0, which 
immediately implies that X2 = 1 since that maximizes I]{f )'E{g) without changing F,{fg). 

To summarize we are down to /([E]) = ^([E]) = 0, /([Sl^y) = g{[S]ai) = gUSja^) = 0, 
/([S](ti) = /([SJctz) = aii^las) = 1> and /([A]) = ^([A]) = 1, so for positive correlation we 
need 

a > {a + + SCT2)(a + s^-J, 

which is true if ae > {s^^ +50-2)- The right hand side is clearly maximized when (T3 = 1 
(since by assumption si is the maximum of si, 52, and S3), and we obtain ae > si(s2 + S3) 
as claimed. □ 

Remark (a). It is clear that the standard FKG technology does not extend to the present 

case. Indeed, if we view our system as {0, 1}"^, but restrict our attention to measures which 
assigns weight zero to the double edge configurations, then it is obvious that the FKG 
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lattice condition fails for any such measure. On a slightly more subtle level, regarding 

{[A], [S]i, [S]2, [SJa, [E]} as simply a partially ordered set with lattice structure given by 
X\/Y = inf{Z\X ^ Z and y ^ Z} and X A y = sup{Z\Z ^ X and Z ^ Y}, it is not 
hard to see that the FKG lattice condition holds whenever ae > siS2- This is in apparent 
contradiction with (the necessity part of) Theorem (2.4). However, the connection between 
the lattice condition and positive correlation hinges on the fact that the lattice satisfies 
distributivity, which is a property that our lattice lacks, as [S]i = [S]i A ([S]2 V [SJs) / 
([S]iA[S]2)V([S]iA[S]3) = [E]. 

Remark (b). We observe that ae > 2s^ implies that a*e* > 2(s*)^ by (2.5), (2.6), and 
(2.7), so the *-map takes the region of positive correlation into itself. 

Remark (c). It is noted that for independent bonds, at density p, the condition ae > 2s^ 
is well-satisfied. But supposing we write 

e={l-pf{l-t), (2.9) ^2^ 

s=p{l-pf{l-t), (2.10) [32b" 

and 

a= (/ + 3p2(l-p))(l-t)+t (2.11) [32^ 

(as we will have occasion to do when we discuss magnetic systems) and again consider, 
with the obvious interpretations, our old eight configurations. Then it is clear that the 
correlations between bonds are positive if and only if i > 0. However, our condition ae > 2s^ 
is satisfied for values of t which are considerably negative. 

2.3 Phase Diagram 

PD Theorem 2.5. Consider the correlated percolation model on the triangular lattice as defined 
previously which has parameters e, s and a; the parameters are described by points in the 
ae-plane. Supposing that ae > 2s^ , then in the region 

2^/2 

a + e > ro = 



3 + 2^2' 
the following hold 

(1) The region a> e has a ( unique ) infinite cluster. 

(2) The region a < e has no infinite cluster and is characterized by exponential decay of 
correlations. 

(3) The line a = e has no infinite cluster of either type but power law (lower bounds) on 
the decay of correlations. 

These results are summarized in Figure 1. 
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Figure 1: Phase diagram for percolation problem on the triangular lattice; variable s suppressed. The line 

a = e is the self dual line. The curve ae = 2s^ separates the regions with and without positive correlations. 
Within the region of positive correlations, a > e is the percolation phase, a < e non-percolating with 
exponential decay of connectivities and percolation of the dual model. These phases are divided by the 
self dual line, where there is no percolation of either type and critical behavior is observed. Some of these 
results may be extended out of the region of positive correlations by monotonicity. 



Proof (sketch): We will be brief since the major ingredients are transcriptions with minimal 

modifications of the well-known results from standard percolation theory. Our setup will be 
as follows: we will fix the value of a + e, denoting this by r, and write a = Xr, e = (1 — A)r, 
< A < 1. We will denote by Ac(r) the (purported) threshold above which there is 
percolation (Notwithstanding, we do not "yet" know that there will be percolation even 
if A = 1). Notice by Proposition (2.1) and its corollary that, in these circumstances, the 
duality takes A to 1 — A. 

Our first claim is that the result on the exponential decay of connectivities below thresh- 
old applies whenever r > [24], [25]. The starting point is the adaptation of Russo's formula 
[29] to the current situation. For an increasing event A, a triangle t is pivotal if, when empty, 
the event A docs not occur but if fully occupied then it does. Denoting by Pj-.x the prob- 
ability measure with parameters a = Xr, etc, and E^^a for the corresponding expectation, 
the modification of Russo's formula is easily shown to be 

= rBrA\SA\), 

where \SA\ denotes the number of pivotal triangles for the event A. 

Next, we denote by An the event that the origin is connected by occupied bonds to the 
boundary of a "ball" of radius n. It is clear that the basic "chain of sausages" picture holds 
in this context (with paths of bonds replaced by clusters of triangles) and at the endpoint 
of each sausage, a pivotal triangle. We note that for the present setup two events are said 
to occTir disjointly if they are determined on the configurations in disjoint sets of triangles. 
Thus, using the more general Reimer's inequality [27] in place of the van den Berg-Kesten 
inequality one can follow the standard derivations to obtain 

Er,x{\SA„\) > - const. 

Thereafter, some tedious analysis shows that if at some Aq, P.,.,Ao(^n) ^ then for all 
A < Ao, 3^* > such that Fr,\iAn) < e~"^('*'''')"; in particular there is exponential decay of 
connectivities. However, standard 2D arguments show that once the direct model has rapid 
decay of correlations, the dual model percolates. (E.g, if there is no connection between 
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points on the x-axis with x < —L and points with x > +L than some dual point with 
x-coordinate in the vicinity of the gap is connected to infinity). 

Using duality this immediately implies that Ac < 1/2: Any other possibility would 
imply percolation of the dual model at values of A greater than 1/2 which, by duality, 
implies percolation at A-values less than 1/2, contradicting the possibility of any other 
possibility. 

For general values of parameters, the results of [6] apply which rules out the possibility 
of multiple infinite clusters (of the same type). In the region of positive correlation (r > ro) 
the results of [28] and [18] (see also the proof by Zhang, 1988, unpublished) demonstrates 
that infinite clusters of the opposite type cannot coexist. This implies that there cannot be 
percolation of either type on the self-dual line, i.e. that Ac > 1/2 so that Ac = 1/2. 

Finally, to prove power law lower bounds on the decay of correlations, we observe that 
for appropriate rectangles of length-scale L, there is either a left-right crossing by the 
direct bonds or a up-down crossing by the dual bonds, so that without loss of generality 
the crossing probabilities are of order unity uniformly in L. Standard arguments (see e.g. 
Theorem 2 in [11]) can then be used to demonstrate power law lower bounds. □ 

Remark. Our assumptions of positive correlations and that si = S2 = S3 are the ingre- 
dients needed to use the Zhang (and [28], [18]) arguments. Without these assumptions we 
cannot mathematically rule out the possibility of percolation before or at the self-dual point 
with unique infinite clusters of both types. In the independent case, coexisting clusters were 
ruled out in [16] using direct (Kesten-style) arguments. It is conceivable that these argu- 
ments could be modified to the present case but we make no specific claims. Nevertheless, 
some of the isotropic results can be extended outside the regions of positive correlations by 
domination arguments: 

Corollary 2.6. In the region a > ro/2, e < ro/2 the relevant (percolative) conclusions 
of Theorem, 2.5 hold while in the region a < ro/2, e > ro/2 the relevant non-percolative 
conclusions of Theorem 2.5 hold. 

Proof: Consider a point with parameters a > ro/2, e < ro/2 which is not covered in 
the previous theorem. Such a point can be joined by a horizontal line to a point in the 
percolative region described in Theorem 2.5. For all intents and purposes, the new measure 
is obtained from the known percolative measure by replacing empty triangles with singly 
occupied triangles: Explicitly, the measure in question stochastically dominates a measure 
with the stated percolative properties. The conclusion follows since the two claims about 
the regions a > e, a + e > ro may be phrased in terms of the events: 

(1) The existence (wpl) of an infinite cluster and 

(2) Uniqueness of said cluster. 

The first is manifestly increasing while the second is equivalent to the absence of an infinite 
cluster of the dual type, hence also increasing. The region a < ro/2, e > ro/2 is handled 
similarly. □ 
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3 Generalized Star— Triangle Relations — Random Cluster Measure 
3.1 Graphical Weights and Spin Systems 

We start in this subsection with the random cluster models - a generalization of the usual 
random cluster models which features interactions among certain triples of sites. Here 
we will confine attention only to triples which constitute three vertices of an up pointing 
triangle. 

The random-cluster models are defined by four parameters, e, s, a, and q, and are given 
formally by 

W{ij) a g^HslsHlal^('^)lel^('^)l (3.1) 

where a; is a bond configuration, |s(a;)| denotes the number of triangles with solely one side 
occupied and 13(0^)1 denotes the number of triangles with all three vertices connected, and 
|e(ci;)| the number of empty triangles. It may be assumed, without loss of generality, that 
a+3s+e = 1. Of course as usual the above only makes sense in finite volumes with particular 
boundary conditions; infinite volume problems are extracted via limits. However, as far as 
we are concerned, boundary conditions only enter through the definition of c{uj); once we 
establish the basic monotonicity properties of the model, there are natural dominations in 
both volume and the various parameters s, a and e. Then, the passage to infinite volume 
follows the exactly the same lines as for the usual random-cluster model. Indeed, as far as 
these general matters are concerned we refer the reader to [17] (see also [8] and [7]) where 
the issues have been discussed in some detail. 

It is also clear (see the above mentioned citations) that for integer q greater than one, 
this random-cluster model is the graphical representation of a (formal) Potts Hamiltonian 
with two and three site interactions: 

-pn = ^ J {Sa^ay + Saya. + Sa^a. ) + KSa^aya. , (3-2) 

<x,y,z> 

where the sum is over all generic up-pointing triangles. We assume that J is positive but 
there is, as of yet, no restriction concerning the parameter k. 

For completeness, a quick derivation proceeds as follows: Let A denote a finite collection 
of (up-pointing) triangles and TC\ the restriction of ?^ to A with free boundary conditions, 
and Z\ the corresponding partition function. Then, 

= E ^'"'^^ = E n (^^-^-y + + + 1)(1 + gSa.aya,), 

ca cta <x,y,z>eA 

where S = e"^ — 1 and g = e'^ — 1, and again with no stipulation about the sign of g. 
Multiplying everything out, we get 

CTA <x,y,z>€A 
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where A = 35^ + g{l + S)^ , which we now stipulate to be positive. Notice that we have 
dehberately failed to distinguish terms corresponding to products of two, versus three Kro- 
neckcr deltas. Opening up the product and identifying graphical terms in the usual fashion 
we perform the trace to obtain 



(3.3) partition 



where uj denotes a bond configuration restricted to five possibilities on each triangle as 
described in the previous section. Since everything is positive, the objects in the above 
summand represent weights for the configurations lo. For convenience, we can multiply the 
above by an overall (irrelevant) factor and then, by suitable redefinitions of parameters, we 
have our weights in the form of (3.1). 

We remark that the more standard decomposition into eight configurations per triangle 
would, as can be checked, lead to positive correlations if and only if g > 0. Indeed, g/{l + g) 
corresponds exactly to the parameter t which was discussed in equations (2.9)-(2.11). As 
we will show in Theorem 3.4 below, the present system provides a great deal more leeway. 

Remark. Finally, it is worth a reminder that as far as the spin systems are concerned, 
most quantity's relevance can be read directly from the graphical problem [1], [13]. In 
particular (at least in the realm of positive correlations), percolation is synonymous with 
magnetization, while the absence of percolation implies unicity among the possible limiting 
Gibbs states. 



3.2 Duality Relations and Self-Dual Curve 

Theorem 3.1. For the random cluster measure as defined in the previous section, the 
duality relations are given by 

e* V a / 

and 



a 

e* 



2 

q^e 



The self dual curve, obtained in the above by setting a = a* , e = e* and s = s* is then 



a = qe. 

Remark. We note that the above corresponds exactly to equation (15) in [35]. 

Proof (sketch): To derive the duality relations, we make use of Euler's formula, which, 
as usual, has to be interpreted in the context of specific boundary conditions. And here we 
have the additional step of integrating out the middle spin to return to the triangular lattice. 
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However, with careful consideration of the situation at the boundary, dual measures may 

be identified in finite vohinic. Specifically, if A consists of nothing more than N connected 
up-pointing triangles with free boundary conditions, then the dual model will consist of 
the corresponding down-pointing triangles with fully wired boundary conditions. Other 
scenarios at the boundary can be treated in a similar fashion; we will be content to proceed 
formally. But before we begin there is yet another technical difficulty: Our three-body 
interactions do not distinguish between triangles with two or three edges occupied; in order 
to use Euler's formula we must take this into account, so we set the convention that all 
three-body interactions have all three edges occupied. Now, finally, we have: 

Thus if a;* is the standard dual (on the hexagonal lattice) we have: 

where |e(a;*)| corresponds to the number of empty triads, etc. Finally, integrating out all 
middle spins, we obtain: 



W{lo) DC q' 



o,q<-') (^y^^^ 



\<^*)\ 



eq 
a 



Ja(c.*)| 



|aK)! 



qmvar 



Here we have used the fact that the empty configuration on the triad has four connected 
components while that on the triangle when the middle vertex is integrated out has only 
three, so we must compensate a factor of q for each e(a;*), yielding the q~^a. The weights are 
now in the form of equation (3.3). Derivation of the self-dual curve is now straightforward. 

□ 

Simple algebra now gives: 
Corollary 3.2. For A > 0, the regions ae > Xs^ are invariant under the *-map. 



3.3 Positive Correlation 

Our proof of positive correlations will concern N triangles with configurations of the type 
described and measures determined by the weights given in (3.1). For the purposes of this 
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measure 



ive correlation 



proof, we make no restrictions on the geometry of the triangles: they need not represent a 
subset of the triangular (or any other planar) lattice. In general, sites can belong to any 
number of triangles, but if a pair of sites belong to two distinctive triangles, the associated 
bonds can appear twice. In addition, we will need to consider different sorts of boundary 
conditions on our N triangles; these will, generically, be denoted by F. These T conditions 
are the identification of sets of points which are considered to be "preconnected" (even if 
no bonds are present). In particular, the specification of T provides us with a precise notion 
of c{uj) and, for all intents and purposes, determines the geometry of the collection. 

Definition 3.3. Fix a, e, s with a + 3s + e = 1. Let A be a fixed set of vertices in the 
triangular lattice corresponding to M triangles which we label ti,t2, ■ ■ ■ ,tM- Let Tm = 
{[E], [S]i, [S]2, [SJs, [A]}^ denote the set of configurations on {ti,t2, ■ ■ ■ ,tM}- Let F denote 
an arbitrary wiring on A, then for lo G Tm, 



where c{uj, F) now denotes the number of connected components determined by the wiring 
F as well as the configuration uj. Now for N < M, and uj G Tjv, we let 

M—N times 

denote the measure on those N triangles obtained from the corresponding weight. 

Remark. The main thing to remember from the above definition is that we are working 
with some a priori A and all the vertices of A are taken into account in the term c(a;,F); 
this will become important later in the section when the structure of the weights actually 
come into play. Needless to say, we will be interested (for the purposes of induction) in an 
N which may be envisioned as far smaller than M; indeed, for finite M there is no difficulty 
with the immediate passage M — oo. 

Theorem 3.4. Let /j^ denote the measures as described above with q > 1 and ae > 2s^. 
Then for all N and all wiring boundary conditions T, these measures have positive correla- 
tions. 

The idea is to proceed by indTiction on the number of triangles which we regard as 
embedded in the larger space of M triangles, — M of which are automatically empty. 
We will need the strong inductive hypothesis that has positive correlations for all 

possible wirings F. For the case N = 1, there are clearly only five possible outside wirings: 
no vertices are connected, the vertices corresponding to side one (respctively two and three) 
are connected, and finally all three vertices are connected; we denote these wirings by S, 



(3.4) weights 
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Si S2, S3, and A, respectively. The all wired case, namely /j.'f is exactly the case proved in 
Theorem 2.4. Let us quickly dispense with another example, /xf ^ . Here we see 

fif'iiSl) = zs. 



Meanwhile, for A; = 2, 3, 
and finally 

/xf([A]) = ^^, Mf([E]) = ^e, 

where we use the notation [A], [S]i, ... ,[E] to denote the relevant corresponding events and z 
is a normalization constant. The necessary inequality /xf^( [A] )/Ltf^([E]) > /xf^([S]i)(/Ltf'([S]2)+ 
/xf'([S]3) follows readily from ae > provided q > I. The other cases are just as easily 
demonstrated and we may consider the base case to be established. 

We make use of two key ideas in the forthcoming inductive proof. The first is a general- 
ized version of the lattice condition. Indeed, whenever the underlying space is the product of 
linearly ordered spaces, the lattice condition is entirely equivalent to the minimalist version: 

Kr?,a,6')-K^,a',6')' ^ ' ^ 

where the a's and 6's represent the configuration at any two coordinates, rj is all other 
coordinates and a > a' and b > b'. Crucial to our argument is that despite the absence (or 
inapplicability) of the full lattice condition, an analogue of (3.5) nevertheless holds. The 
second key idea is a slight generalization of Proposition 2.22 in [22] which is the statement 
that a convex combination of two measures with positive correlations itself has positive 
correlations if one of the measures FKG dominates the other. We state and prove these as 
our next two lemmas below. 

1P2T Lemma 3.5. Let /x^ be defined as above with q > 1. Then an analogue of (3.5) holds 
for jji^, provided the separate increases pertain to different triangles. E.g., «/Tiv-2 is the 
configuration on the first N—2 triangles, and we have T^-i, T'^_^, Tjv, T'j^j as configurations 
on the last two triangles with Tjv-i ^ ^jv-i '^''^'^ '^n ^ T'j^, then 

Proof: Examining the ratios in the statement above, a quick glance at (3.4) reveals that 
all the "prefactors" cancel on both sides of the purported inequality, leaving only the q- 
dependent terms. Since q > 1, the above amounts to a special case of 

C{uj \/ri)+ C{uj Ari)> C{uj) + 0(1]), 

which has been proved in complete generality in numerous places (e.g. [1]). □ 
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Ill 



increasing 



Lemma 3.6. Let {H, >z^) and (/C, >z^) be finite partially ordered sets. Let ji he a probability 
measure on H and for each ri £ Ti, let v-ri be a probability measure on K. It is supposed that 
fi has positive correlations, that for each rj, the measure has positive correlations and 



furthermore, if rji y rj2, then Vj^^ > 



Then 



has positive correlations. 

Proof: Let / and g be increasing functions on /C. Then 



to 



= ^fi{v)B,^{fg) 

r)£H 

>^/z(r?)E,^(/)E,^(5). 

It is observed from the hypothesis that Ej^^ (/) and Ej^^ (5) are increasing in ry and the result 
fohows from the positive correlation of ^u. □ 

Now let us informally proceed with an inductive proof. In what is to follow we assume 
that/ and g are increasing functions on N triangles, Tjv-i always denotes the configuration 
on the first N — 1 triangles and Tjv € {[A], . . . , [E]} a generic state of the N^^ triangle. We 
condition on the state of the last triangle, and according to Bayes' formula, we get 

Tn 

where MatiAjv restriction of fi^ to the last triangle. 

As far as the first — 1 triangles are concerned, we can apply the inductive hypothesis 
to conclude that the measures |Tjv) has positive correlations, since the conditioning, 

along with F, give us some wiring scenario for these triangles. So (appealing to Lemma 
(3.6)) we will be done if we can show that (i) E^(/|Tjv) and E^(5|r^) are increasing in 
Tn (i.e. ij^{—\Tn) > )U^(— |Tj^) whenever Tjv >: Tj^), and (ii) the measure /i^^iA has 

FKG ' 

positive correlation. These are the topics of yet the next two lemmas. 
Lemma 3.7. Let f and Tjv be as described and define 

Ft, = EUf\TN). 



Then Ftj^ is an increasing function. 
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Proof: Suppose Tjv >r TJ^. Then we note that, as in the standard argument, Lemma (3.5) 
impUes that 

(^(Tjv-i) = -jr— — p— — (3.6) g 

is an increasing function of Tjv-i = (Ti, ...,Tjv_i). We aim to show that 

(/|Tjv) > (/it;,). 



We have 



Tjv-i 

/iiv(rjv) 

since / is increasing and T/v > T'j^. Now the last expression can be rewritten as 
E /(T^-i.r^) ^^^^r';:^Vf^V (T^-i) = E^(/0|T^), 



Tjv-i 



which by induction is greater than or equal to E^{f\T'p^)'E^{(j)\T'^). Thus, concatenating 
the above expressions, we have 



Er(/|T^)>Er(/|T^)Er(0|T; 



N) 



= ( 2^ /(Tjv-i,rjv) — — )( 2^ — - — ) 

sinceVx ^^(^ . ^-^-f^) = 1. □ 

poscor Lemma 3.8. Let ^^|^^ denote the measure jj,^ as described above restricted to the N^'^ 
triangle. Then ^^|/^^ has positive correlation. 

Proof: We will again make use of Lemma 3.6, so we write 

B 

where B represents the total wiring conditions outside the N^^ triangle due to the initial 
wiring condition F and the outside configurations, Tjv-i- However, the overall effect of F 
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and Tat-i is to produce one of the five types of wiring on a single triangle - a situation with 
which we are familiar - and henceforth wc may assume B G {A,Si,S2,S3,£}. 

We note that for each B, ^Ajv positive correlations (the subject of the base case). 
Next, we claim that B ^ B' implies that fiff > fif^. This follows from the observation 

FKG 

that less wiring on the outside produces more factors of for the weights of the higher 
configurations (see the Remark after Definition (3.3)). Explicitly, it can be checked that for 

where D (which depends on B' and B) is a decreasing function of T^. Thus we have verified 

two of the three hypotheses of Lemma (3.6). 

We are down to the last hypothesis; here we will need to write fi^{B) in a more explicit 
form. Note that by induction ii^_-^^{B) has positive correlations, so we seek some relationship 
between fJ.^{B) and fjJ^_i{B). To do this we must exploit the "almost" product structure of 
the weights (3.4) from which our measures came. So wc first let Z]v-i(S) denote the weight 
of observing B, before the introduction of the A''*'^ triangle, and let Zjj_-^^ = Xli? ■^a^-i('^) 
denote the overall normalization factor, so that \b = Zn-i{B)/Z'^_^ = ii^_^{B). Next 
we may write Zn^i{B) = I]t;v-i lirjv-iur=BWj(TAr-i), where W^(Tiv-i) is the weight of 
observing the configuration Tjv_i as given by (3.4). Similarly, if Zjv(S) denotes the weight of 
observing B given the iV*'* triangle, then Z^iB) = Etat-i lT;v-iur=s (J2tn (Tjv-i, Tat)^ 
Comparing the previous two expressions and referring back to Definition (3.3), it is not dif- 
ficult to see that 

ZN{B) = n^ZN-i{B), 

where (up to a factors of e) = (^ + y + e), n^. = (| + s + y + e), and = (a + 3s + e) 

- which happens to be one. Thus, letting Zjf = J2b ^n{B), we arrive at 

It is noted that the factor Zjf_-^^/Zj^ is independent of the wiring B, T^, etc. Thus by 
Theorem (2.4) all we need to show is that {n^XjO{n^X£) exceeds n_j^ A^^ (ra^^ '^''^'Ss^Sa) 

- or whatever ordering combination maximizes the latter object. To this end, let cr be a 
permutation on three letters such that A^^^ > A^^^ and Xs„^ > Xs„^ ■ Our last hypothesis 
will be verified if we can show that 

{ngX£){n^XA) > {ris^^Xs^Jin^^^^Xs^^ +"-5,3^5.3)- 

To this end, we first observe that the induction hypothesis implies A^Ag > A^^^^j ('^5r(2) + 
-^5^(3)) for any permutation on three letters r: On general grounds this is true because of 
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the similarity between the outside wiring space and the inside configuration space. But, 

to proceed formally, let / and g be the increasing functions of the outside wiring such 
that fi£) = g{£) = 0, f{A) = g{A) = 1, = 1 - <7(5,(i)) = 1 and = 

1 — = 0, i = 2, 3. Then by the fact that hn-i has positive correlation, we indeed get 

^A^£ ^ ''^<St(i) (^•St(2) +-^'Sr(3) )• On the basis of this inequality we only need that n^n^ > 
(since n, = n,, = n„ ), i.e. we need that, 

/ ~ \ / CL 3(5 > , (X 25 > o 

(a + 3s + e)(^ + — + e)>(- + s + — + e)^ 
q q q 

Now if one multiplies and compares terms, one has an expression which involves {q — 1) 
times a quantity which is "easily positive" . 

We have verified all three hypotheses of Lemma (3.6) and can therefore conclude that 
/ij^l^^ has positive correlation. □ 

Proof of Theorem (3.4): As already remarked, we will (again) use Lemma 3.6. Explicitly, 
we apply Lemma (3.6) with 7i = {[A], [S]i, . . . , [E]} (corresponding to configurations on the 
A'"*'* triangle) and ^ = fi^, , and /C = {[A], [S]i, . . . , [E]}^ (corresponding to configurations 

on all N triangles) and = /i^(— |t/). The three hypotheses of the Lemma are verified by 
the induction hypothesis and Lemmas 3.7 and 3.8. □ 

We conclude this section with a brief discussion on infinite-volume limits: In the region 
of positive correlations, more wiring leads to a higher measure. Thus, for free boundary 
conditions (the restrictions of) measures increase with increase volume and for fully wired 
boundary conditions, they decrease. So, for a nested sequence of volumes which exhaust 
the lattice, well-defined infinite- volume limits - which are independent of sequence - exist. 
Furthermore, as mentioned earlier, wired and free measures may be dually identified in 
finite volume. Thus, in turn, we may identify the dual of the infinite volume free measure 
as the wired measure and vice versa. 



3.4 Phase Transitions 

In this subsection, we establish results on phase transitions in the g— state Potts/random 
cluster models under consideration. Here, unlike in the percolation case, we cannot establish 
with certainty whether the transition is continuous or discontinuous. Moreover, for the 
continuous cases, our statements will be considerably weaker than our Theorem 2.5 since 
much of the technical artillery (e.g. the van den Berg-Kcstcn inequalities) do not apply. In 
particular, in the continuous case, we cannot even prove that the percolation/magnetization 
transition actually occurs on the self-dual line. Nevertheless, critical behavior is established 
for self-dual points which are also points of continuity, the subject of our first proposition: 
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Proposition 3.9. Consider the random cluster model on the triangular lattice as defined 

by (3.1) and satisfying ae > 2s^. Then at any self-dual point a = qe which is a point of 
continuity of the bond density the following hold: (1) The percolation probability vanishes 
and (2) there are power law lower bounds on the correlation functions. 

Proof: Much of the proof can be transcribed directly from our proof of Theorem 2.5 and 
as for the rest, similar arguments have appeared before ([3], [11]), so we will be succinct. 
The first statement follows from the results in [19] which, under the conditions of positive 
correlations and 2D symmetries, forbids coexisting infinite clusters of the opposite types. 
Thus, in any realization, there is either no percolation of either type or there arc separate 
states (depending on how the infinite-volume limit was constructed) which have and don't 
have percolation. However, this latter happenstance, by appeal to Strassen's Theorem 
[30] implies that the distinctive states have different bond densities which would imply 
a discontinuity in the bond density. For the second statement, routine arguments which 
may be traced back to [1] imply that the limiting random cluster measure is unique and 
therefore may be identified with the dual measure; on this basis the rest of the argument 
follows mutatis mutantis from the proof of Theorem 2.5 for percolation (again see [11]). □ 

Finally we show that in the region of positive correlation, any discontinuity in bond 
density must occur on the self-dual curve: 

Proposition 3.10. Consider the random cluster model on the triangular lattice as defined 
by (3.1) and satisfying ae > 2s^. Then any discontinuity in the bond density must occur on 
the self-dual curve as given by a = qe. 

Proof: Our proof is a variation of the one found in [3] ; here we unfortunately do not have 
a convenient family of curves which are nicely preserved under the duality relations. We 
will work with A and S parameters given in (3.3); suppose at (^O)S'o) - with Aq > 2Sq 
- we have a discontinuity in the bond density. Let A > 1 and consider the curve C = 
{(2AS'q,5'o) : a > 1}. We note that along C the measure is FKG increasing with increasing 
A. Next let A'^^, A^ and A^ denote the corresponding values of A at which the curve C 
intersects the self-dual curve, the percolation threshold, and the discontinuity, respectively. 
We aim to show that A"^-^ = = A^. 

Let Ci denote the part of C which is below the self-dual curve and similarly let Cu 
denote the part of C which is above the self-dual curve. Since it is not the case that 
{CiY = Cu, we need to define two new curves to work with: Let Cy = Ci U (Q)* and 
Cp = C„ U (Cu)* , and we parametrize Cy by Xy and Cp by Ap with the requirement that 
A = Ay on Ci and A = Ap on C„ (and extending in the obvious fashion). We remark that 
with these parametrizations, Cy and Cp are FKG increasing in A^ and Ap by duality (or by 
explicit computation): E.g., on Cy PI Ci the measures are clearly FKG increasing; on the 
other hand, this implies the measures corresponding to the dual parameters - these lie on 
{CiY and are parametrized by A^ > A^^ - are decreasing in Xy for Xy < Xy^, and hence 
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increasing m Aj, for Xy > X^^. Now if we let Xy, Xy^ and A^ denote the corresponding 
values of Aj^ at which the curve C intersects the self-dual curve, the percolation threshold 
and (should it exist) the discontinuity, respectively. Similarly we define A^, A^^ and A^ 
for Cp. Then A^^ = A^^ = A^-^. 

First we show that A^ > A*^^: If this is not the case, then for A'^^ > A > A^ the 
direct model is percolating in the wired state. Note that this A corresponds to a A^ in 
our new parametrization. At the dual value, A*, we would then have dual percolation in 
the state with free boundary conditions. However, the dual model in the wired state at 
parameter A "dominates" the dual model in the free state at parameter A*, and hence there 
is dual and direct percolation at A (e.g. in the wired state), which is a contradiction of 
[19]. Next we can easily show that A^ < A^: This is because a discontinuity in the bond 
density implies non-uniqueness of the limiting measure and hence, ultimately, percolation. 
Finally, we must have A^ < A'^^: Towards a contradiction assume that X^ > A^-°; this 
implies that A^ actually exists and is equal to A^. Next note that the same argument that 
showed A^ < A^ also shows A^ < A^. Since we have a discontinuity in the direct model 
if and only if we have a discontinuity in the dual model, we have another discontinuity at 
A* < X^^ < Xp, a contradiction. □ 

4 Conclusion 

We have described a Potts/random cluster model on the triangular lattice with three-body 
correlations. By introducing a reduced state space, the duality relations are easily derived. 
It is noted, in the context of spin systems, that the purely ferromagnetic region of parameters 
is not mapped into itself under dTiality. More generally, in the q > 1 random cluster 
models, when we consider the full state space, the region which has positive correlations is 
not mapped into itself. However, for the reduced case, necessary and sufficient conditions 
for positive correlations are derived which are invariant under duality and include a larger 
portion of the original parameter space. Under the conditions of positive correlations, 
for percolation and for values of q where there are discontinuities, it is proved that the 
transition occurs at the self-dual point; if there is no discontinuity, self-dual points admit 
critical behavior. On the basis of exact solutions [4] it has been argued that the dividing line 
is g = 4, similar to the situation on the square lattice. The advantage of the current random 
cluster formulation is that this hypothesis can be tested numerically using cluster methods; 
e.g., the algorithms in [31], [23] and [10] can be readily adapted. While we have no reason 
to doubt the results in [4] in this case, for a related model with three-body interactions on 
the square lattice, there is some evidence pointing to the reduction of the dividing q. In 
any case, although we will not discuss details, it should at least be possible to prove that 
for large q there is a discontinuous transition. Here certain modifications will be needed to 
adapt the methods of reflection positivity to the present case, which may very well be the 
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subject of a later paper. 
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